第 39 卷 第 6 期 应 用 A 学 学 报 ,Yol39No.6 
2022 年 12 月 Chinese Journal of Applied Mechanics ~ ^ NE 1 FA Dec. 2022 


文章 编号 :1000- 4939 ( 2022) 06-1065- 11 
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HO 要 :无 网 格 形 函 数 的 非 多 项 式 特 性 导致 梯度 计算 复杂 耗 时 ,同时 高 斯 积分 方法 不 满足 积分 约束 
条 件 , 因 此 传统 伽 辽 人 金 无 网 格 法 难以 达到 理论 收敛 率 。 再 生 光 滑梯 度 的 构造 特点 使 其 自然 满足 积 
分 约束 条 件 ,并 有 效 地 避免 计算 无 网 格 形 函 数 的 直接 梯度 ,因而 具有 高 效 和 精确 的 特点 。 为 了 探 完 
再 生 光 涓 梯度 无 网 格 法 的 动力 特性 ,本 研究 构造 了 基于 再 生 光 滑梯 度 理论 的 伽 辽 金 无 网 格 法 的 动 
力 分 析 方 法 ,详细 研究 了 再 生 光 滑梯 度 无 网 格 法 的 动力 计算 精度 ,包括 频 散 特性 、 自 由 振动 和 时 程 
动力 分 析 。 再 生 光 滑梯 度 无 网 格 法 采用 再 生 光 滑梯 度 替代 传统 的 无 网 格 形 函 数 梯度 ,由 于 其 本 身 
与 积分 约束 条 件 的 内 在 一 致 性 ,直接 采用 基 函 数 对 应 阶 次 的 低 阶 高 斯 积分 方法 对 质量 和 刚度 矩阵 
进行 数值 积分 , 即 可 保证 最 优 收敛 率 和 精度 。 理 论 分 析 与 数值 计算 结果 均 表明 ,再 生 光 滑梯 度 无 网 
i 格 法 的 频 散 特性 、 频 率 收 伊 率 和 时 程 动力 计算 精度 ,都 明显 优 于 采用 高 阶 高 斯 积分 方法 的 传统 无 网 
格 法 。 
关键 词 :无 网 格 法 ;动力 分 析 ; 再 生 光滑 梯度 ;收敛 性 ;精度 
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Characteristics of dynamic meshfree analysis with reproducing 
kernel smoothed gradients 
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Abstract: The gradient valuation of non-polynomial type meshfree shape functions is usually very complex 
and costly , which leads to an undesirable fact,i. e. , even very high order Gauss quadrature rules cannot 
ensure an optimal convergence of Galerkin meshfree formulation. The reproducing kernel smoothed gradi- 
ents of meshfree shape functions naturally meet the Galerkin integration constraint , which completely avoid 
the time consuming computation of meshfree shape function gradients , and thus is highly efficient and ac- 
curate. In order to assess the dynamic performance of reproducing kernel smoothed gradient meshfree algo- 


rithm ,this work presents a dynamic meshfree analysis with reproducing kernel smoothed gradients and the 
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convergence and accuracy are particularly investigated in detail, including the dispersion analysis , free vi- 


bration analysis and transient analysis. Since the present algorithm inherits the integration consistency of 


Galerkin formalism, the relatively low order Gauss quadrature rules corresponding to the basis order in 


meshfree approximation are capable of producing the optimal convergence with superior accuracy. Both 


theoretical and computational results consistently demonstrate that, regarding the dispersion analysis, fre- 


quency convergence as well as transient response, the reproducing kernel smoothed gradient Galerkin 


meshfree formulation gains much more favorable solution accuracy for dynamic meshfree analysis, com- 


pared with the conventional meshfree scheme using high order Gauss integration. 


Key words :meshfree method ; dynamic analysis ; reproducing kernel smoothed gradient ; convergence ; accu- 
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无 网 格 法 '" 将 计算 区 域 离散 为 一 系列 节点 , 仅 
利用 节点 的 位 置信 息 即 可 构造 高 阶 光滑 ,全 域 协调 
的 形 函 数 , 因 而 被 广泛 地 应 用 于 大 变形 分 析 :9 、 薄 板 
者 高 阶 问题 计算 "各 向 异性 材料 热 弹 性 问题 分 
WESS \ 结 构 振 动 分 析 """ 等 。 虽然 伽 辽 金 无 网 格 法 
共有 和 角 棒 性 好 、 精 度 高 的 特点 ,但 是 一 般 情 况 下 ,无 
网 格 形 函 数 为 有 理 式 ,并 且 背景 积 分 单元 与 形 函 数 
影响 域 无 规则 相交 ,因而 实际 计算 中 为 了 得 到 理想 
的 请 算 结果 ,每 个 背景 积分 单元 需要 采用 高 阶 高 斯 
积 宗 。 另 一 方面 ,高 阶 高 斯 积分 方法 的 采用 又 会 严 
恒隆 低 无 网 格 法 的 计算 效率 "* ,阻碍 了 无 网 格 法 
模 琢 大 规模 问题 。 因 此 ,发 展 高 效 精确 的 伽 辽 金 无 
网 格 数 值 积分 方法 便 成 为 无 网 格 法 研究 的 一 个 重要 
方面 。 针 对 上 述 问题 ,许多 学 者 提出 了 各 具 特 色 的 
DERIK. PU, Chen 等 55 系统 提出 了 伽 辽 
金光 网 格 法 的 线性 积分 约束 条 件 ,发 展 了 高 效 的 且 
能 够 适用 于 线性 基 范 数 的 稳定 节点 积分 方案 。 针 对 
TARDI R, Liu 等 "9 研究 了 一 种 径 向 点 插值 积 
分 方法 ,但 该 方法 的 计算 结果 依赖 参数 的 选择 。 
Duan 等 "1 拓展 了 线性 积分 约束 条 件 ,提出 了 一 致 积 
分 法 ,该 方法 满足 了 高 阶 积分 约束 条 件 ,但 其 需要 在 
每 个 背景 积分 单元 求解 方程 组 来 构造 形 函 数 梯度 。 
Chen 等 59] 在 彼得 罗 夫 伽 辽 金 框架 下 ,通过 重 构 特 
定 的 权 函 数 ,建立 了 能 够 满足 任意 阶 积 分 约束 条 件 
的 数值 积分 方法 ,但 是 彼得 罗 夫 伽 辽 金 列 式 会 产生 
非 对 称 刚 度 和 矩阵。Hillman 等 "对 稳定 和 变 分 一 至 
的 节点 积分 方案 进行 了 系统 的 研究 ,并 将 其 应 用 至 
冲击 破坏 问题 。 

Wang 等 "基于 应 变 光 滑 方法 和 理 查 德 外 推 
法 ,发 展 了 一 种 具有 二 阶 准确 特性 的 拒 套 子 域 积 4 
高 效 无 网 格 法 。 随 后 , Wang 等 站 进一步 建立 了 具 
有 显 式 积分 特性 的 , 且 适 用 于 任意 次 基 函 数 的 再 生 


光滑 梯度 无 网 格 法 理论 框架 ,简洁 起 见 , 也 称 其 为 再 
生 梯 度 法 。 该 理论 框架 满足 伽 辽 金 无 网 格 法 的 积 
约束 条 件 与 无 网 格 形 函 数 的 梯度 再 生 条 件 , 具 有 高 
精度 ,高 效率 的 特点 ,并 能 够 保证 伽 辽 金 无 网 格 法 的 
一 致 收 化 性 ”1。 在 此 基础 上 , 匡 俊 超 等 站 基于 赫 
林 格 - 赖 斯 纳 ( Hellinger-Reissner) 变 分 原理 , 发展 了 
一 种 具有 变 分 一 致 特性 的 高 效 无 网 格 强 制 边 界 条 件 
施加 方法 。 该 方法 保证 了 离散 系统 的 稳定 性 ,同时 
不 依赖 任何 人 工 参 数 。 另 一 方面 ,无 网 格 法 与 机 器 
学 习 的 结合 近期 也 得 到 了 关注 。 例 如 , 陈 健 等 .1 基 
于 循环 卷 积 神经 网 络 提出 了 一 种 无 网 格 动力 分 析 代 
理 模 型 ,其 不 仅 可 以 加 速 动力 分 析 的 计算 过 程 ,同时 
可 以 保证 数值 计算 精度 。 针 对 热 参 数 辨识 问题 , 梁 
钰 等 “| 采用 径 向 基 函 数 训练 热 参 数 ,构造 了 一 种 参 
数 空间 与 时 域 空间 双向 压缩 代理 模型 ,具有 能 够 快 
速 搜索 任意 参数 的 时 域 响应 特点 。 

值得 指出 的 是 ,再 生 梯度 法 "虽然 实现 了 精确 
和 高 效 的 静 力 分 析 , 但 仍 缺 乏 对 其 动力 特性 的 研究 。 
针对 该 问题 ,本 研究 分 析 了 再 生 梯度 法 的 动力 特性 ， 
包括 频 散 特性 .自由 振动 和 时 程 分 析 。 结 果 表 明 ,与 
采用 传统 形 函 数 梯度 的 伽 辽 金 无 网 格 法 相 比 ,再 生 
梯度 法 不 仅 具 有 更 优 的 频 散 精度 和 收敛 特性 ,而 且 
所 需 高 斯 积分 点 数 较 少 ,计算 效率 更 高 。 


1 无 网 格 近 似 与 再 生 光 滑梯 度 理论 


1.1 再 生 核 无 网 格 形 函 数 


无 网 格 法 通过 在 分 析 区 域 2 和 边界 厂 上 布置 
一 系列 节点 {x,21 来 实现 空间 离散 ,其 中 NP 代表 
无 网 格 节点 数量 。 每 个 无 网 格 节 点 x 对 应 的 形 函 
数 为 W, (x) ,影响 域 为 supp (x,)。 考 虑 任意 变量 


投稿 网 站 :http://cjam. xjtu. edu. en 微 信 公 众 号 :应 用 力学 学 报 


第 6 期 


“其 对 应 的 无 网 格 近 似 函数 w(z ,人 为 
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Hop ic WAHE d, DIWA s i 的 系数 。 
根据 再 生 核 近似 理论 ”” ,无 网 格 形 函 
为 如 下 形式 。 
W(x) =p (x,-x)c(x) eg, (x, -x) (2) 
其 中 :c(x) 为 待定 系数 向 量 ;p(x) 为 p 阶 单项 式 基 
向 量 , 且 
p(x) 911 uy, a y o Y Ositjsp(3) 
其 中 wp,(Cz -x) 为 附属 于 节点 x 的 核 函 数 ,其 影响 
域 的 大 小 影响 域 尺 寸 s, 确定 ， 核 函数 及 其 影响 域 
v, (x) -(p.(x, -x)A 


u'(x,t) = 


数 可 假设 


> 由 式 (4) 和 (6) 的 表达 式 可 以 看 出 ,无 网 格 形 函 
类 3 般 为 有 理 式 ,因而 其 梯度 计算 复杂 耗 时 。 为 了 
避免 直接 计算 无 网 格 形 函数 梯度 ,提高 计算 效率 , 文 
献 E21] 提 出 了 适用 于 任意 次 无 网 格 基 函 数 的 再 生 
浴 滑 梯度 构造 理论 。 


m CE PUAN 


NN 假设 场 变 量 u(x,t) 为 任意 的 p 阶 多 项 式 , 则 其 
is, (x,t) 可 以 表示 为 
: u(x,t) =a; (1)q(x) (7) 
Pia ,为 任意 系数 向 量 ; ;4(X) 为 (p -1) 阶 的 单项 
式 花 向 量 , 即 q(x) =p "(x)。 由 文献 [15,17] 可 
知 (p -1) 阶 积分 约束 条 件 为 
© [Pan = [ vandr- f Yaa (8) 
由 式 (8) 可 见 , 对 于 了 次 基 函 数 , 伽 辽 金 弱 形 式 
所 采用 的 数值 积分 方法 只 有 精确 满足 (P — 1) 阶 积 
分 约束 条 件 , 无 网 格 数值 解 才能 重 现 对 应 的 多 项 式 
精确 解 。 
根据 再 生 光 滑梯 度 理论 ”  ， 
水 数 类 似 , JE P R JÉ PK 
可 表示 为 如 下 形式 。 


与 再 生 核 无 网 格 形 
3, 的 再 生 光 滑梯 度 更 ，， 


P(x) 24 (x)c (x) 9 (x) (9) 


其 中 :ec, 为 待定 系数 向 量 ; p (x) 为 核 函数 ,这 里 
取 为 
1, xef 


$0 «|, A^ (10) 
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!(x)g(x, -x) *p'(x, -x)A;! 
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的 大 小 共同 决定 了 无 网 格 形 函数 的 局 部 紧 支 性 和 光 
滑 性 。 方 便 起 见 ,本 研究 均 采 用 3 次 B FEZR EZ ERE 
进行 计算 ,二 维 问题 的 影响 域 为 张 量 积 形 式 和 矩形 域 。 
将 一 致 性 条 件 引 入 式 (2) , 即 可 得 到 未 知 系数 
向 量 e (x) ,进而 可 得 到 无 网 格 形 函 数 表 达 式 为 
P(x) =p (0)A  (x)p(x, -x)g(x, -x) (4) 


AF A (x) NIE EIERE, H. 
A(x) = 2 p(x, -x)p'(x, - x)e(x, - x) 
(5) 
对 式 (4) 进行 微分 , 即 可 得 到 无 网 格 形 函数 的 


一 阶梯 度 为 
(x)e(x, =x) *p'(x, -x)A (x)p (x -x))p' (0) 
(6) 
其 中 :0Q2, HERES H UX Qe =N 的 积分 单元 ; 
Ne 表示 积分 单元 的 总 个 数 。 图 1 给 出 了 再 生 光 滑 
梯度 无 网 格 法 采用 的 三 角形 背景 积分 单元 。 不 失 一 
般 性 ,这 里 以 三 次 基 薄 数 为 例 详细 阐明 再 生 光 滑梯 
度 构 造 过 程 。 当 采用 三 次 基 添 数 p(x) 时 ,gqg (x) 的 
q(x)2(L,w,y,v y uy) xeQ; (11) 
将 式 (11) 代 入 积分 约束 条 件 式 (8) ,可 得 


| WadQ- g1,1212, NP (12) 
Qc 


&j-| andr -| Pga (13) 
Te Qc 


再 用 式 (9) 定 义 的 再 生 光 滑梯 度 多 ,替换 到 式 (12) 


中 的 标准 梯度 V, ,, 可 得 ce; = Gi'g 5, 其 中 G6 为 再 生 
光滑 梯度 的 矩 量 矩 阵 , 且 
G, = fa q'dQ (14) 


便于 后 续 数 值 积 分 ,如 图 1 所 示 , 这 里 将 三 角形 
积分 域 投 影 至 参数 空间 。 在 投影 后 的 三 角形 积分 域 


内 采用 高 斯 积分 即 可 求解 得 到 g ? ,具体 表达 式 为 


Nc 


z Y > V,(xcs)q(xog) nJ (Xon) wen 一 


CTI cB ACT o) 
Nc 
V,(xq) q i(Xa)J(GXa) Oi 
CTI GL ZACOG) 
(15) 
其 中 :A 为 背景 积分 单元 内 部 或 边界 上 的 高 斯 积 4 
点 总 数量 ;xcen .wc8 为 背景 积分 单元 边界 上 的 高 斯 积 
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分 点 位 置 和 相应 的 权重 iX .wa 为 背景 积分 单元 内 
的 高 斯 积分 点 位 置 与 配套 的 权重 ;J 为 背景 积分 单 
元 联系 物理 和 参数 空间 的 雅克 比 和 矩阵 行列 式 。 注 意 
到 式 (14) 中 的 G6 可 以 直接 解析 计算 ,无 需 采用 数值 
BUT, 

结合 式 (9) 及 (15 ) , 即 可 得 背景 积分 域内 的 再 
生 光 滑梯 度 表达 式 为 

W(x) = (x) G'g i, xee (16) 

Hrs 16) 的 再 生 光 滑梯 度 表 达 式 可 以 看 出 ,再 
生 光 滑梯 度 通过 直接 构造 得 到 ,避免 了 标准 无 网 格 
梯度 的 复杂 计算 , 显著 提高 了 梯度 计算 效率 。 图 2 
给 出 了 一 维和 二 维 情况 下 采用 三 次 基 函 数 构造 的 无 
网 格 形 函 数 标准 梯度 与 再 生 光 滑梯 度 。 


f ; 一 — 
Y 


S 图 1 


三 角形 背景 积分 单元 示意 图 


c Fig.1 Schematic diagram of triangular 


Q background integration cells 

需要 指出 的 是 ,上 述 再 生 光 滑梯 度 仅 用 于 刚度 
矩阵 的 构建 ,而 对 于 后 处 理 过 程 中 的 应 变 和 应 力 计 
算 , 仍 采用 式 (4) 和 (6) 的 无 网 格 形 函 数 及 其 梯度 。 
这 样 做 的 主要 原因 是 再 生 光 滑梯 度 并 非 全 域 连续 也 
数 , 若 采 用 其 进行 后 处 理 还 需要 进一步 采用 其 他 的 
全 域 应 力 光滑 方法 。 因 此 ,一 种 比较 简便 的 后 处 理 
方法 是 基于 再 生 光 滑梯 度 所 求 得 的 高 精度 节点 系 
数 ,采用 全 域 光 滑 的 标准 形 函 数 梯度 直接 计算 应 变 
和 应 力 结果 。 

此 外 ,注意 到 无 网 格 形 函 数 再 生 光 滑梯 度 的 影 
响 域 与 标准 梯度 的 影响 域 不 完全 相同 。 这 是 由 于 再 
生 光 滑梯 度 是 基于 背景 积分 单元 构造 的 ,与 无 网 格 
形 函 数 影响 域 相交 的 背景 积分 单元 组 成 了 再 生 光 滑 
梯度 的 影响 域 ,所 以 再 生 光 滑梯 度 的 影响 域 一 般 略 


ChinaX 


Vv 合作 期 刊 第 39 着 


大 于 无 网 格 形 函数 标准 梯度 的 影响 域 。 当 采用 线性 
基 汪 数 时 , 式 (16) 即 退化 为 稳定 节点 积分 采用 的 光 
滑梯 度 ™ 。 


2310 zo 
00 02 04 06 08 LO 00 02 04 0.66 08 10 
x x 


| p-3 


Tw 
y 0000 € 


习 2 一 维和 二 维 三 次 基 无 网 格 形 函数 的 
标准 梯度 与 再 生 光 滑梯 度 


Fig.2 Standard and reproducing kernel smoothed 


gradients for one-dimensional and two-dimensional 


meshfree shape function with cubic basis 


2 再 生 光 滑梯 度 无 网 格 法 的 离散 方程 
及 频 散 精度 


2.1 离散 方程 
不 失 一 般 性 ,考虑 如 下 典型 波动 方程 , 即 
u(x,t) =c [u(x,t) Ft Cet) ] +b(x,t), xeQ 


u(x,t) =g(x,t) " xel* 


en. Vu(x,t)=t(x,t), xel" 
u(x,0) =u (x), xeQ 
u(x,0) =u (x), xeQ 


(17) 
其 中 :e 为 波 速 ;Y 为 空间 梯度 算 子 ;2 DU UIN; IU 和 
分 别 代表 自然 和 强制 边界 ;1 和 e 为 对 应 的 给 定 
边界 值 ;n = {n,n 为 边界 的 外 法 向 ; 场 变 量 w 上 
部 的 点 和 下 标的 逗号 分 别 表示 时 间 和 空间 导数 。 
对 于 式 (17) 所 定义 的 波动 方程 ,其 等 效 积 分 弱 
形式 为 
| suin +e | Vôu - VudQ = | 8ubdQ + | Suar 
NQ N Q r 


(18) 
其 中 6 为 变 分 符号 。 将 式 (1) 代 入 式 (18) 中 , 即 可 
得 到 无 网 格 运动 离散 方程 为 
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第 6 其 
Md «Kd - f (19) 
其 中 ,MK 和 了 分 别 表示 一 致 质量 和 矩阵、 刚度 矩阵 
和 力 向 量 ,其 具体 表达 式 为 
M, = | V,(x) V,(x)dQ (20) 


K; = e| CQ)" Vw,x)dQ (21) 


pe I:ZOLUE NES idr (22) 
TURN (LT) PIS V (x) 替换 为 再 生 光 滑梯 度 
V V. (x) , 即 可 得 到 再 生 梯 度 法 的 刚度 矩阵 表达 
式 为 
Ky =ef (VÉ A(x)) VE (x)d2 (23) 
> 与 式 (15 ) 相同 , 3X E Nil BER RE E — $ p gt B E 


的 于 算 均 采 用 高 ! 斯 积分 方法 。 因 此 , 式 (20) 和 (23 ) 
NAEMIAUTON 


Nc 


Z7 B 2: b» P xa) Pia) J Xc) 9 
x 4 C-1GI-A(QQ) 
£ J Nc pe as — 
KON N= > V V,(xq) V V (Xa) J (Xc) ® c 
= C-1GI ZA(OQ) 
(25) 
n: A AERE AE 


^ [f ,无 网 格 形 函 数 的 非 插值 特点 导致 其 无 法 
像 有 限 元 法 一 样 直接 施加 强制 边界 , 因此 需要 辅助 
其 它 方 法 来 进行 施加 强制 边界 ,例如 拉 格 朗 日 乘 子 
活 C 吧 函数 法 或 尼 效 ( Nitsche ) 方法 等 “2 。 本 研究 
利用 再 生 光 滑梯 度 的 内 在 特点 ,采用 文献 [24] 提出 
的 赫 林 格 - 赖 斯 纳 变 分 一 致 方法 来 施加 边界 条 件 , 但 
这 里 给 出 了 该 方法 更 广义 的 加 权 残 值 构造 格式 。 根 
据 该 方法 , 场 变 量 ”和 其 梯度 ,为 两 个 独立 变量 ， 
此 时 考虑 引入 边界 条 件 的 等 效 积分 弱 形式 为 


0 = | u(cu,-*b-u)dQ- 
0 ' 


Í ôu cn, (u - g)dr - 
Te 
| Sulu in; - dr (26) 
Tt 
对 式 (26) 引 入 分 部 积分 ,有 
0 =f 8uc^u ,n,d - | ôu cu dQ + | 6ubd0 - 
re 0 0 


| 5viaQ = f uc naar + 


付 赛 赛 ， 等 : 再生 光滑 梯度 无 网 格 法 动力 特性 研究 ChinaXiV 合 1 期刊。 1069 


f Bu c^ n,gdT s]: SuidT (27) 
进一步 利用 自然 边界 D 和 强制 边界 D* 的 互补 
关系 , 即 T= 了 TAT', 有 


| 8uc'u nidT = | 5ucu nd - | uc u inidT 
Te r n 
= | ôu cu dQ + Í 8uc'u dQ - 
o 0 
Í 8uc^u nidT (28) 
Tt 
Í ôu ic niudT = Í ôu ĉnudr - | ôu ic nud 
re r n 
- | Bu sc udQ + | ôu ,c'u AQ - 
0 0 
Í ôu ^ nud (29) 
I' 


将 式 (28 ) 和 (29) 代 入 式 (27) ,并 利用 6u FI ôu ， 
的 独立 性 ,可 得 如 下 的 等 效 积 分 弱 形 式 , 即 
| Vôu - VudQ = 


Í (n - vou) udT = (V> V&u)udQ + 
I' n 
Í (n * Vôu) gdr (30) 
TE 
| uian + | buc (n ， vu)dr - | Suc ( V- vu) dn 
Q rt 


= Í ôubdN + | SuidT (31) 
Q I'" 


E epic 只 分 单元 Qe 中 ,将 无 网 格 场 
离散 式 (1) 以 及 梯度 离散 式 (7) 代 入 式 (30)， 


ag 


(gS -g5)d, + Gz g? (32) 


其 中 De H Qe 的 边界 ; g EIR (13) gE lg? 见 式 


(33) 
gf yea = | 
利用 式 (32) , HEE u 的 离散 梯度 可 以 表示 为 

u' -2q'a - FY, 


qn,edT (33) 


= Y, 4, T q” Ge gi 
(34) 


其 中 更 ， 兄 式 (16) V, =q" GZ gi. 
将 式 (1) 和 (34) 代 入 式 (31) , 即 得 考虑 边界 条 
件 的 系统 离散 方程 "为 
Md +Kd=f (35) 


K=KR+K+kK,f=F+f+f (36) 
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sth M AR R24) ~ (25); F 65k (22), K K 
和 力 向 量 了、f 中 的 元 素 见 式 (37)。 


| Y, n'è V Wr - | (V w)'env,r 
k , 


-| ( V V,)'cngdl 
" 


(37) 


Ky =| (FY) "envar 
LES (38) 


lox | - (VV,)'cngdr 


HopVw,-|v,.. T d ETIL BERAI 


DE K 5 f imh m c9 PRG SU E PEU RO E 
Ris 与 尼 效 法 六 不 同 , 这 里 的 稳定 项 不 需要 人 工 
IA. 

加 对 于 时 间 域 离散 ， 本 研究 采用 HHT ( Hilbert- 
Hüuahg transform) 77 dz 79! ,该 方法 仍然 采用 Newmark 
淡 的 预测 -矫正 算法 ,位移 和 速度 系数 向 量 预测 和 矫 
下 阶段 的 表达 式 分 别 为 

AJI) 预测 阶段 


2 
c — 


. (39) 
Z va mv, * CAD(1 -y)a 
^2) MEE EE 
c =d  * BAD? 
一 5 N ELANA (40) 
Van —-Vyn 


+1 +yAt Ql 


n :7?=d,xa=y;86 和 7 是 积分 参数 ;Al: 是 从 tn 到 
tw 时 刻 的 时 间 步 长 ,下 标 “m” 表示 相应 时 刻 的 量 。 

根据 HHT 771 7^ ,时 间 - 空 间 全 离散 的 运动 方 
程 为 


Ma,,;*(1*o)Kd,,,-oKd,-f,,, (41) 

fusa =f nsa) tuus tuu tat (42) 

其 中 wa 为 HHT 方法 的 参数 ,当选 取 wel -1/3,0], 

B-(1-20)/4,y 2 (1 -2a)72 时 ,该 时 间 积分 方法 

条 件 稳定 。 将 式 (40) 代 入 式 (41), 即 可 得 到 加 速 
度 系数 向 量 的 更 新 方程 为 

[M+ (1 +a)B(At) Ka,, 
=f, +akd, - (1*o)Kd,,, 
由 式 (43 ) 解 得 a ,1 后 , 即 可 通过 式 (40) 计 算 ,i 时 
刻 的 速度 和 位 移 系 数 向 量 。 


(43) 
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2.2 频 散 精度 分 析 


频 散 特性 是 度量 离散 动力 系统 精度 的 一 种 有 效 
方法 。 在 频 散 特性 中 ,研究 对 象 为 自由 振动 ， e. 
外 力 和 边界 条 件 的 影响 ,这 时 运动 方程 式 (35 ) 只 
2 200 00 RSS 


和 再 生 光 滑梯 度 组 成 的 刚度 矩阵 玉 ,对 应 的 自由 振 
动 方程 为 
- (w^) 2, Mud; * Y Kyd,-0 


Ip A3 
其 中 :Sj = (JIsupp(x;) msupp(zy) Ø}; w" 表示 频 
率 数 值 解 。 
对 节点 系数 d, (1) 引入 简 谐 波 假定 ,有 
d,(t) =doexp[i(kxj ky, -wt)] (45) 
式 中 :du 为 波幅 ;c= V -lk k, 分 别 为 x.y 方向 上 
的 波 数 。 将 式 (45) 代 入 式 (44) ,可 得 相 速 度 e, DN 


(44) 


(kh) p K y expli (Ex, *tkx)] 


JeSyj 


w 
w > Mh” "expl t (E s, + kx)) ] 


JeSy 


(46) 


式 中 :wo c V+ 名 为 真实 频率 ;h 为 节点 间距 。 另 
= REGI LIBEROS 


e= Cap -1=% -1 = 00") 


式 中 7 表示 频率 的 收敛 阶 次 。 

图 3 给 出 了 二 维 问 题 的 频 散 计算 结果 ,其 中 
"GEI" 表示 采用 7 个 高 斯 积分 点 的 传统 高 斯 积分 无 
网 格 法 ，RKCS 表示 再 生 光 滑梯 度 无 网 格 法 。 表 1 
列 出 了 各 种 方法 在 每 个 背景 积分 单元 内 所 采用 的 高 
斯 积分 点 数目 。 图 3 的 频 散 分 析 结 果 表 明 ,再 生 梯 
度 法 明显 优 于 传统 无 网 格 法 的 精度 。 


(47) 


90 一 GL3 
PEE =O 
— RKGS 


180 


210 330 


270 
图 3 不 同方 法 的 频 散 精度 对 比 


Fig.3 Comparison of dispersion analysis results 


by various meshfree formulations 
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3 数值 算 例 


本 节 通 过 两 个 典型 算 例 ,对 基于 一 致 质量 矩阵 
的 再 生 梯 度 法 的 动力 特性 及 计算 精度 进行 系统 分 
析 。 算 例 中 传统 无 网 格 法 与 再 生 梯度 法 所 采用 的 高 
斯 积分 点 数 见 表 1 。 
表 1 不 同方 法 使 用 的 高 斯 积分 点 数量 


Tab.1 The number of Gauss points used by various 


meshfree formulations 


计算 JE pat 
模型 阶 次 


标准 高 斯 
积分 点 数 


高 阶 高 斯 RKGS 
积分 点 数 积分 点 数 


1 


ea -1.0 -0.8 -0.6 -04 -0.2 


-1.0 -0.8 -0.6 -0.4 -0.2 
lgh 
图 4 一 维 杆 问题 一 、 二 阶 频率 收敛 结果 对 比 


Fig.4 Comparison of convergence of the first and second frequencies for the 1D rod problem 


为 了 更 好 地 验证 再 生 梯 度 法 的 动力 特性 , 进 一 
步 分 析 一 维 杆 问题 的 时 程 响 应 。 计 算 中 按照 精确 


X 


f(x) e sin( F )sin( 到) 来 设 定 初 位 移 和 初速 
BE ,无 网 格 节点 数 为 11,HHT 时 间 积 分 方法 的 参数 


ESI 
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3.1 一 维 杆 问 题 


考虑 一 维 两 端 固定 杆 的 自由 振动 问题 。 不 失 一 
般 性 , 取 杆 长 二 =5, 横 截面 面积 4、 材料 密度 p 及 杨 
REE E 均 为 1。 对 于 该 问题 , 其 频率 解析 解 为 
w; =imc/L, 其 中 波 速 。= VE/p。 无 网 格 计算 中 采用 
11 ,21 41 个 均 布 节点 进行 收 伍 率 分 析 , 其 中 二 次 和 
三 次 基因 数 的 影响 域 大 小 分 别 为 2.5 和 3.8。 图 4 
给 出 了 标准 伽 辽 金 无 网 格 法 和 再 生 梯 度 法 分 别 采 用 
二 次 与 三 次 基 函 数 时 前 两 阶 频率 的 收敛 结果 。 结 果 
表明 ,再 生 梯 度 法 的 频率 计算 结果 对 于 二 次 和 三 次 
基 函 数 分 别 可 以 达到 4 次 和 6 KAREAR ， 
而 标准 伽 辽 金 无 网 格 法 即使 采用 高 阶 高 斯 积分 仍然 
无 法 达到 相应 的 最 优 收敛 阶 次。 同时 ,再 生 梯度 法 
在 采用 较 少 高 斯 积分 点 情况 下 ,其 频率 计算 精度 也 
远 优 于 标准 伽 辽 金 无 网 格 法 ,充分 说 明了 再 生 梯 度 
法 在 频率 计算 方面 的 精确 性 和 高 效 性 。 


-1.0 -0.8 -0.6 -04 =02 


-1.0 -0.8 -0.6 -04 -02 
lgh 


a = -0.3, 时 间 积 分 步 长 Ar =0. 0125, R5 给 出 了 
FHKE KZA 3 种 方法 分 别 在 x =1.0 和 
x=2.5 位 置 的 位 移 时 程 曲线 误差 , 可见, 再 生 梯 度 法 
的 动力 时 程 分 析 误 差 明 显 低 于 采用 低 阶 ( GI-2 ) 和 高 
阶 高 斯 积分 (G1-3 ) 无 网 格 法 的 动力 时 程 分 析 误 差 。 
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图 5 一 维 杆 问 题 动力 分 析 位 移 误差 对 比 


Fig.5 Comparison of displacement error for the dynamic analysis of 1D rod problem 


3.2. 二 维 膜 问题 
考虑 二 维 四 边 固定 方 膜 自由 振动 和 动力 分 析 问 


了 均 布 和 非 均 布 无 网 格 离散 情况 下 二 维 膜 问题 的 
一 、 二 阶 和 二 、 三 阶 频率 的 收敛 率 。 可 见 ， 无论 是 
均 布 还 是 非 均 布 无 网 格 离散 ， 再 生 梯 度 法 均 能 达到 


题 总 何 及 材料 参数 为 : 边 长 了 ,厚度 ~\ 密 度 及 面 
内 张力 7 均 取 相应 的 单位 值 1,c = T pr) ,其 频 


最 优 收敛 率 2p， 对 于 二 次 和 三 次 基 消 数 即 为 4 和 
6。 相 较 而 言 ,采用 3 点 和 7 点 (二 次 基 范 数 )、6 点 


看 节点 离散 模型 进行 分 析 ， 有 具体 的 节点 数 依次 为 
121: 441 和 1681。 无 网 格 形 函 数 采 用 二 次 和 三 次 


和 12 点 


LN 


EA 六 E mj x N y 
基调 数 时 的 影响 域 分 别 为 28 813.6, 图 7~8 给 出 无 网 格 法 。 
T INNNNNNNNNNNNNNNNNNAN 
INNNNNNNNNNNNNNNNNNNN 
INNNNNNNNNNNNNNNNNNNN 
INNNNNNNNNNNNNNNNNNNN 
NINNNNNNNNNNNNNNNNNNN 
INNNNNNNNNNNNNNNNNNNN 
NNNNNNNNNNNNNNNNNNNN 
INNNNNNNNNNNNNNNNNNNN 
INNNNNNNNNNNNNNNNNNNN 
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NSI 
N 


441 nodes 


121 nodes 


点 (三 次 基 函 数 ) 高 斯 积分 的 传统 无 网 格 法 呈 


图 6 二 维 腊 问题 的 均 布 和 非 均 布 无 网 格 离散 模型 


现 出 精度 低 、 收 敛 率 低 等 问题 。 


注意 到 对 于 三 次 基 


函数 ,频率 计算 结果 甚至 出 现 不 收敛 问题 。 
一 维 计算 结果 类 似 , 再 生 梯度 法 对 于 二 维 问题 的 频 
率 收敛 率 和 精度 同样 明显 优 于 采用 高 斯 积分 的 传统 


1 681 nodes 


因此 ,与 


Fig.6 Uniform and non-uniform meshfree discretizations for the 2D membrane problem 
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Fig.7 Comparison of convergence of the first and second frequencies for the 2D 
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图 8 二 维 腊 问题 非 均 布 无 网 格 离散 二 三 阶 频率 收敛 结果 对 比 


Fig.8 Comparison of convergence of the second and third frequencies for the 2D membrane problem 


with non-uniform meshfree discretizations 
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对 于 该 膜 问题 ,同样 进行 了 动力 时 程 分 析 , 计算 
中 根据 解析 解 ux Lt) =2sin( T) sin( 7) sin (V2 rt) 


选取 初 位 移 和 初速 度 。HHT 方法 中 的 参数 
a= -0.3, 时 间 积 分 步 长 At =0. 001。 图 9 给 出 了 该 
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SP 


N 
交通 过 再 生 光滑 梯度 理论 构建 的 无 网 格 形 函 数 光 
PE ,满足 积分 约束 条 件 ,消除 了 无 网 格 形 函 数 全 
元 得 时 的 传统 梯度 计算 ,提高 了 无 网 格 法 的 精度 与 
效率 。 本 研究 以 再 生 光 滑梯 度 理论 为 基础 ,详细 讨 
论 字 再 生 光滑 梯度 无 网 格 法 的 动力 离散 方程 ,并 通 
过 芯 林 格 - 束 斯 纳 变 分 方法 施加 强制 边界 条 件 ,保证 
了 窗 散 系统 的 变 分 一 致 性 。 最 后 ,研究 了 再 生 梯 度 
法 的 动力 特性 ,包括 频 散 特性 、 自 由 振动 和 时 程 动力 
At 

结果 表明 ,再 生 梯度 法 的 频 散 精度 明显 高 于 采 
用 高 斯 积分 的 传统 伽 辽 金 无 网 格 法 。 在 计算 精度 方 
面 ,由 于 再 生 梯度 与 积分 约束 条 件 具有 内 在 一 致 性 ， 
因而 采用 较 少 的 高 斯 积分 点 即 可 达到 频率 的 最 优 收 
敛 率 。 对 于 二 维 二 次 基 函 数 ,每 个 背景 积分 单元 仅 
采用 3 点 高 斯 即 可 达到 4 次 收敛 率 ,而 传统 无 网 格 
法 即便 在 每 个 背景 积分 单元 采用 7 个 高 斯 积分 点 ， 
仍然 无 法 达到 最 优 收敛 率 , 日 频 率 计 算 精 度 也 远 低 
于 再 生 光滑 梯度 无 网 格 法 ; 当 采 用 三 次 基 函 数 时 , 采 
用 12 点 高 斯 积分 的 传统 无 网 格 伽 辽 爹 法 甚至 出 现 
频率 不 收敛 的 情况 。 时 程 动力 分 析 的 结果 同样 证 
明 , 再 生 光 滑梯 度 无 网 格 法 的 动力 计算 精度 远 优 于 
传统 伽 辽 金 无 网 格 法 。 
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模型 在 采用 二 次 基 辆 数 和 121 个 均 布 节点 离散 情况 
下 ,点 (0.2,0.5) 和 点 (0.4,0.4) 的 位 移 时 程 曲线 误 
差 对 比 结果 。 动 力 分 析 结 果 再 次 验证 了 再 生 梯度 法 
与 采用 高 斯 积分 传统 无 网 格 法 相 比 的 精度 和 效率 
优势 。 
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图 9 二 维 膜 问题 动力 分 析 位 移 误 差 对 比 


Fig.9 Comparison of displacement error for the dynamic analysis of 2D membrane problem 
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